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Abstract. The degree of linear polarization has recently been found 
to show wide random variations over the solar disk. These variations 
are presumably at least partly due to fluctuations in the flux density of 
turbulent photospheric magnetic fields and associated variations in the 
degree of Hanle depolarization. In order to understand the origin of such 
large scale fluctuations of the turbulent magnetic flux density we develop 
a phenomenological model to calculate the spatial Fourier spectrum of 
the fluctuations of turbulent magnetic fields in the solar photosphere and 
convective zone. It is found that if the model parameters are fitted to 
turbulence closure models and numerical experiments the characteristic 
scale of the fluctuations is by about an order of magnitude larger than 
the turbulence scale (the scale of the granulation), owing to the more 
effective quenching of small-scale fluctuations by turbulent diffusion. 



1. Introduction 

In a highly conductive and strongly turbulent plasma like the solar photosphere 
random small-scale deviations of the flow field from mirror symmetry are ex- 
pected to incessantly generate small magnetic flux loops of random orientation. 
This is a miniature analogue of the large-scale dynamo process occurring in a 
globally non-mirror-symmetric flow, for which reason it is known as small-scale 
dynamo action (Leorat, Pouquet, & Frisch 1981, Kida, Yanase, & Mizushima 



1991, Petrovay & Szakaly 1993). Owing to the random orientation of the loops, 
no net large-scale field will arise, but a non-zero mean magnetic energy den- 
sity B 2 and mean unsigned flux density \B\ results. Most of this turbulent flux 
resides in magnetic structures with scales much smaller than the characteristic 
scale of the turbulent velocity field (which is the granular scale of ~ 1000 km in 
the solar case). Traditional Zeeman magnetography is "blind" to these fields as 
the net circular polarization of the mixed polarity small-scale field cancels out 
in a resolution element. 

It has long been recognized that the Hanle effect (magnetic depolarization 
of linearly polarized radiation) may offer a way to detect the turbulent fields 



In: High Resolution Solar Physics: Theory, Observations, and Techniques (Rimmele 
T., Balasubramaniam K. S. & Radick R. eds.), Proc. 19th NSO/SP Workshop, Conf. 
Series Astron. Soc. Pacif., in press 



1 



Fluctuations of Turbulent Magnetic Field 



2 



(Stenflo 1982| ). The observational study of turbulent fields has however been 
hampered by the shortage of observational data, by the lack of a reliable radia- 
tive transfer theory for polarized radiation in magnetic fields and by the fact 
that beside turbulent fields, the Hanle effect is also due to resolved magnetic 
elements (network, ephemeral active regions) and to the overlying canopy fields 
(especially for lines formed higher in the atmosphere). Nevertheless, in recent 



years important advances have been made in all these areas (Faurobert-Scholl 


|1993l Faurobert-Scholl et al. 1994, Bianda, Solanki, & Stenfloj 


1998|, Bianda, 


Stenflo, & Solanki 


1995, Landi Degl'Innocenti 1996, Trujillo Bueno & Landi 


Degl'Innocenti 1997 


, Landi Degl'Innocenti 1998). 



On of the most intriguing recent discoveries is the finding of Stenflo, Keller 



fc Gandorfer ( 1998 ) that the degree of linear polarization Q/I shows large ampli- 
tude random variations over the solar disk. While, as mentioned above, canopy 
fields, network elements and ephemeral active regions may also contribute to the 
observed Hanle depolarization, it is still likely that at least part of this observed 
spatial variation is due to the presence of similar fluctuations in the flux density 
|J3| of the turbulent photospheric magnetic field. 

The presence of fluctuations should not come as a surprise from a theoretical 
point of view. After all, in the turbulent solar photosphere any physical quantity 
should show fluctuations, and even variations of an amplitude comparable to the 
mean value are commonplace. What is more surprising is the spatial scale of 
the fluctuations. While the observations have a resolution of about 1" along the 
slit, the dominant variations seem to occur on the much larger scale of ~ 10 4 km. 
V ariations on the granular scale are of much smaller amplitude. (See e.g. Fig. 3 
in Stenflo, Keller, fe Gandorfer 1998| .) The real theoretical challenge is therefore 
to understand how the dominant scale of fluctations of turbulent flux density 
can be so much larger than the turbulence scale? 

One popular explanation for the existence of large-scale photospheric struc- 
tures (such as supergranulation) is that they are the "imprints" of processes 
going on deeper down in the convective zone where the characteristic scales (de- 
termined by the pressure scale height H) are larger. Alternatively, it is of course 
also possible that the large scales are due to some local photospheric process 
like an inverse cascade. In order to resolve this problem, we need a model for 
the generation and transport of turbulent magnetic flux that takes into account 
both the generation and saturation processes constituting the small-scale dy- 
namo and the turbulent transport of |-B| throughout the underlying convective 
zone. In the following such a model will be presented. 



2. The spectrum of fluctuations: a linear model 

The evolution equation for the unsigned flux density \B\ of the turbulent field 
should read something like 

(\B\/B y \B\/t- (1) 
s * ' 

nonlinear saturation 
q > 



dt\B\= V((3V\B\) + |£|/t + 

linear generat 
(small-scale dynamo) 



turbulent diffusion linear Keneration 
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As the transport of the highly intermittent magnetic field in a turbulent plasma 
is due to advection of flux tubes irrespective of their polarity, the transport 
terms are expected to be identical to those for the signed field B. For simplicity, 
here we consider an isotropic turbulent diffusion with diffusivity (3 = lv/3 (first 
term on the r.h.s.), I, and v being the correlation length and r.m.s. amplitude 
of the turbulent velocity field. The linear generation term corresponding to a 
small scale dynamo would lead to an exponential growth of a homogeneous \B\ 
field to infinity, were it not for a higher order term leading to the saturation 
of \B\ at a finite value Bq, induced by the curvature force (last term on r.h.s.). 



Numerical simulations and closure calculations (|Kida, Yanase, & Mizushima 


1991 


, Nordlund et al. 


1992, 


Durney, De Young, & Roxburgh! |1993 




De Young 


1980 


) show that Bq is about an order of magnitude lower than the equipartition 



flux density: Bq ~ B eq /10. 



As the coefficients 1/t + and 1/t_ are functionals of the turbulent velocity 
field they are expected to show considerable fluctuations around their mean 
values. Hence, the flux density \B\ determined by the equilibrium of generation 
and nonlinear saturation processes will also fluctuate around its mean value 
\B\ ~ Bq. It is plausible to assume 



1/ T _ = 1/ T+ = 1/r = v/l. (2) 

Furthermore, introducing the notation \B\/Bq = 1 + 6, it greatly simplifies the 
treatment if one assumes 6< 1. (In Section 3 we will consider the problem to 
what extent this simplification affects the results.) This allows us to linearize 
equation (1). In the case when b varies much faster with depth than Bq we 
obtain 

d t b = V(/3Vfo) - qb/r + (1/t+ - 1/r-)' (3) 

Now the second term on the r.h.s. describes the net mean restoring effect of 
generation and saturation terms, tending to reduce the fluctuation b, while the 
last term is the fluctuation generation term or forcing term arising owing to the 
turbulent fluctuations in the coefficients. 

We further introduce the (physically plausible) assumption that the Fourier 
spectrum of this forcing term is dominated by those isotropic modes whose 
vertical phase is such that at the depth where the pressure scale height equals 
the inverse of their horizontal wavenumber k they are maximal: 

(1/t+ - 1/t_)' = - h( k > z ) exp[i(fcr + ky + ut)] (4) 
T k 

(Note that throughout this paper we take the logarithmic pressure InP as in- 
dependent variable in the vertical direction, the depth z being just a shorthand 
notation for a function z(lnP), determined by a convection zone model, x and 
y are the horizontal coordinates.) For the solution of equation (3) we take a 
similar Ansatz: 

b = 2_. b(k, z) exp[i(kx + ky + ojt)] (5) 
k 

Let us consider one mode only. (Note that in this case b can always be 
considered real as by virtue of our assumption about the vertical phase of modes 
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b(k,z) may be written as b\{k) cosflnP — lnPo(A;)], and the initial phase in b\ 
can be chosen freely by a displacement of the time axis.) We substitute (4) and 
(5) into (3), simplify, and take the real part. To simplify the notation, from this 
point onwards we omit the hats. With this we arrive at 

-d z ((id z b) + {2k 2 (3 + q/T)b = bf/T (6) 

This equation determines the Fourier amplitude b(z, k) of each mode of hori- 
zontal wavenumber k in the spectrum of turbulent magnetic field fluctuations. 
Vertical diffusion is now separated in the first term; horizontal diffusion and the 
restoring force constitute the second and third terms on the l.h.s. 

The r.h.s. is the Fourier amplitude of the forcing, bf is the fluctuation 
amplitude produced in time r by the forcing if other terms were not present. As 
the forcing is due to the action of the fluctuating flow field on the existing 
magnetic field, it is plausible to assume bf/r oc B^/B^u, where B^ is the 
corresponding Fourier amplitude in the spectrum of \B\, and Tk(k, z) is the eddy 
turnover time. This spectrum may be determined from closure calculations, 
numerical simulations and observations. Its simplest representation is by two 
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Figure 1. Vertical profiles of some Fourier modes in the spectrum of 
fluctuations of the turbulent magnetic field. The modes are labelled by 
their value of 1/k in Megameter units. Dashed: Bq 

power laws joining in a peak at k ~ 10/L Thus, we represent the r.h.s. by 

b f (z, k)/r(z, k) = 10 2 / 3 " p (k/kOf, k ~ 1/H(z), (7) 

where p = pi for k < 10/cq and P2 otherwise. (Note that in fact a spectrum with 
two breaks might be more realistic as the spectrum of r is peaked at a lower 
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wavenumber of k ~ l/l; our purpose here, however, is just to present a simple 
example calculation.) On the basis of the observed properties of photospheric 
magnetic fields and motion, p\ ~ 1.5 seems to be a realistic choice. The value 
of P2 will be found to be irrelevant to the solution below. 

For the solution of equation (6) the parameters (3 and r are interpolated 
from a convective zone model, q is evaluated by fitting a solution of equation (1), 
with the diffusive term and the perturbations of the coefficients neglected, to 
profiles of |i?| 2 (£) resulting from the closure model of De Young] ( |198C| ); this yields 
q ~ 0.1. The boundary conditions may be chosen as closed (F = j3 d z b = 0) or 
open {F = vb); experimenting has shown that this does no exert a very strong 
influence on the resulting mode profiles. For the numerical solution of equation 
(6) we take In P as independent variable and use a relaxational method. 

The vertical profiles of some Fourier modes in the spectrum of fluctuations 
of the turbulent magnetic field are shown in Figure 1. As expected, the modes 
with larger horizontal scales have their maxima in deeper layers. 



1 .2 




Figure 2. Fourier spectra of the fluctuations of the turbulent mag- 
netic flux density as functions of the horizontal wavenumber. Solid: 
pi = 1.5, q = 0.1; dashed: p\ = 0.5, q = 0.1; dash-dotted: p\ = 1.5, 
q = 1.0; dotted: p\ = 0.5, q = 0.1 with vertical transport switched off 

Computing a large number of modes and plotting their amplitudes near the 
surface against k yields the Fourier spectrum of the fluctuations in the turbu- 
lent magnetic field (Fig. 2). Experimenting with the parameters in equation (6) 
we find that the choice of p2 is irrelevant for the spectrum, B>o/B eq determines 
its amplitude, while p\ and q determine the shape of the spectrum. A striking 
feature of the spectra in Figure 2 is that their maxima fall to significantly lower 
wavenumbers than k§ = l/l, i.e. the characteristic scale of the fluctuations in 
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turbulent magnetic flux density is much larger than the typical scale of turbu- 
lent motions. The physical background of this phenomenon is that the high 
wavenumber components are more efficiently suppressed by diffusion (2k 2 (5b 
term in eq. (6)). The wavenumber of the maximum increases with q (dashed vs. 
dash-dotted curves), as a stronger nonlinear saturation reduces the role of the 
diffusive terms in the equation. On the other hand, the spectral amplitude at 
even lower wavenumbers depends strongly on the spectral index p\ of the forcing 
(dashed vs. solid curves) — clearly, the shallower the forcing spectrum, the more 
energy is input at the larger scales. 

3. The effect of nonlinearity 

The results presented in the previous section were computed under the assump- 
tion that the fluctuations of \B\ have a small relative amplitude. This assump- 
tion is rather dubious in the light of the large Fourier amplitudes found in the 
model (Fig. 1). In order to have an idea about the extent to which nonlinear 
effects may modify the linear results, in this section we present an alternative 
model that calculates the fluctuating field without the assumption of linearity, 
at the cost of a strong simplification of the geometry: only k = modes are 
considered. This obviously implies that spatial spectra or correlation lengths 
cannot be studied; instead, we will compare the temporal autocorrelation of the 
fluctuations in the linear and nonlinear cases. 
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Figure 3. Snapshots of the solution of equation (7) with q = 1 at 
different instants of time. Crosses: initial profile. 

For k = modes d x = d y = so we write equation (1) in the form 

dt\B\ = d z ((3d z \B\) + \B\/t[1 - (\B\/B y]+ B f /r (8) 

where the last term corresponds to fluctuation forcing by to the turbulent fluc- 
tuations of the velocity field. We model this term as a Gaussian stationary 
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random process with correlation time t(z), vertical correlation "length" 1 (in 
InP units), and a mean displacement of 0.3Bq over r . 

Equation (8) is then integrated numerically starting from an arbitrary per- 
turbed initial state. An example solution is presented in Figure 3. Figure 4 
shows the autocorrelation of the fluctuating flux density as a function of the 
time shift for different cases. It is apparent (dotted vs. dashed curves) that 
the correlation time is primarily determined by the value of the q nonlinearity 
parameter, lower q values corresponding to longer correlation times. This result 
is a close analogue to the findings of Section 2 with respect to the spatial cor- 
relations. On the other hand, switching off the diffusive term in equation (8) 
(dash-dotted vs. dashed curves) has no significant effect, showing that diffusive 
quenching of small spatial scales is not the mechanism responsable for the ex- 
tended correlations here. Replacing equation (8) with its linear equivalent (solid 
vs. dashed curves) does not lead to great modifications in the results. This may 
reassure us to some extent of the reliability of the findings of Section 2 above. 



1 .0 



g °- 8 

O 

$ 

CD 

g 0.6 

e 

o 

"8 

"S 0.4 
s~ 
s~ 
o 
o 
o 

"3 0.2 

0.0 

100 1000 10000 

Time [s] 

Figure 4. Autocorrelation of the fluctuating flux density as a function 
of the time difference at a fixed point in the convective zone in the 
nonlinear model with q = 0.1. Dashed: reference model; dotted: q = 1; 
dash-dotted: reference model with diffusive term switched off; solid: 
linear equivalent 



4. Conclusion 

The results of Section 2 now enable us to answer the question posed at the end 
of the Introduction. Turbulent transport processes indeed result in a dominant 
scale for the fluctuations of the turbulent magnetic flux density that is an order 
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of magnitude larger than the turbulence scale. The physical mechanism behind 
this phenomenon is that small-scale fluctuations are more efficiently damped by 
diffusion, shifting the spectral peak to lower wavenumbers. The low value of the 
q nonlinearity parameter suggested by closure calculations and simulations also 
favors larger dominant scales. 

On the other hand, the possibility that deeper structures are "imprinted" 
on the photospheric pattern can apparently be discarded. Switching off the first 
term (vertical diffusion) in equation (6) does not lead to any reduction of the 
dominant scales (dotted vs. dashed curves in Fig. 2; dash-dotted vs. dashed in 
Fig. 4). 

All this shows that the observations of strongly varying Hanle depolarization 
over the solar disk may be understood from a theoretical point of view even 
if all the depolarization is attributed to turbulent fields. Further advances on 
both the observational and theoretical side may make more detailed comparisons 
between the models and the observations possible, thereby offering the prospect 
of a direct observational diagnostics of the properties of the turbulent dynamo. 
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